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Q^ 1 Abstract 

<i> : 

, We consider the problem of a degenerate electron gas in the background of a 

uniformly distributed positive charge, ensuring overall neutrality of the system, in 
the presence of non-commutativity. In contrast to previous calculations [3] that 
did not include twisted statistics, we find corrections to the ground state energy 
already at first order in perturbation theory when the twisted statistics is taken 
into account. These corrections arise since the interaction energy is sensitive to two 
particle correlations, which are modified for twisted anti-commutation relations. 
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u 1 Introduction 
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The study of non-commutative (NC) geometry and its implications have gained consid- 
erable importance in recent times as these studies are motivated from string theory and 
certain condensed matter systems like the quantum Hall effect. Despite this, the physical 
consequences of non-commutativity remains unclear. In free space (no boundaries) and 
in the absence of interactions, non-commutativity seems to have no observable physical 
consequences [H [2] . On the other hand, in [3] it was found that non-commutativity does 
effect the ground state energy of a degenerate electron gas to second order in perturba- 
tion theory. However, in this computation the role of twisted statistics was not taken into 
account. As is well known by now it is necessary to twist the anti-commutation relations 
in order to restore the Galilean or, more generally, Poincare invariance [5j [6] . Here we 
revisit this calculation taking also due care of the twisted statistics. One expects that 
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the twisted statistics will also have an effect on the ground state energy as it changes the 
two particle correlations [4] and it is well known that the two particle interaction energy 
is sensitive to the two particle correlations. Thus there are two possible ways in which 
non-commutativity may have a physical effect. The first is due to the non-commutative 
nature of space per se and the second the modification of two particle correlations due to 
the twisted statistics required to restore Poincare invariance. 

Here we demonstrate this in the specific setting of [3] , namely, a non-relativistic degen- 
erate electron gas in the presence of a uniform neutralising background charge, interacting 
through a screened coulomb potential. One of the motivations for studying this particular 
system is that this is the typical setting encountered in astro physical or quantum Hall 
systems. 

It is generally believed that the non- commutative parameter 9^ u is of the order of 
the Planck area. It, therefore, seems very unlikely that any experimental signature of 
its presence can be detected by any terrestial observation, unless one finds a way of 
amplifying it by several orders of magnitude. One of the motivations for this paper is to 
study whether this is at all feasible by considering a non-relativistic degenerate electron 
gas, which presumably can contain few orders higher than Avogadro's number (~ 10 23 ) 
of particles, so that the effects of non-commutativity may be amplified. 

The plan of the paper is as follows. We first introduce the basic notations and con- 
ventions of NC geometry in Section [21 In Section [3] we provide a brief review on the 
origin of twisted (anti)-commutation relations. In Section H] we consider the effect of the 
non-commutativity by computing the energy-shift, arising from the screened inter-particle 
Coulomb potential, to first order in perturbation theory. Here we use the technology de- 
veloped by Fetter and Walecka [8j for the same commutative problem [9 = 0). Finally we 
conclude in Section [51 



2 NC geometry 

To fix our notation and conventions, we briefly recall the essentials of a NC geometry. The 
canonical NC structure is given by the following operator valued space-time coordinates, 

[xZ p ,x" op } = i9>» (1) 

Instead of working with functions of these operator-valued coordinates, one can alterna- 
tively work with functions of c-numbered coordinates provided one composes the functions 
(taken to belong to the Schwartz Class of functions i.e. whose derivatives of all order, 
both in configuration and momentum space, vanishes asymtotically) through the Moyal 
star (*) product defined as[7j 

a * /3(x) = [aexpC- 0"" d v )0\(x) (2) 
r = ei?, x = (x°,x\....,x d ). (3) 

With this, the Moyal bracket [x^, x u ]* = (x^ * x v — x u * x^) = iQ^ v is isomorphic to the 
corresponding commutator (pp) involving x£ p . In this paper, we consider only the problem 
of spatial non-commutativity and therefore set 9° l = 0. The Poincare group V or the 
diffeomorphism group D, which acts on the NC space-time R d+l , defines a natural action 
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on smooth functions a G C°°(R d+1 ) as 

(goc)(x) = a(^ _1 x), (4) 
for g £ P or G D. However, in general 

(ga)* e (gP)^g(a*eP), (5) 

showing that the action of the group V or T> is not an automorphism of the Algebra 
A e (R d+1 ), unless one considers the translational sub-group. The Poincare symmetry 
can, however, be restored by the twisted implementation of the Lorentz group, as has 
been shown recently in the literature (5j E]. This in turn implies that the permutation 
group, used for defining bosons or fermions, gets twisted as well in order to maintain 
the statistics operator as a super-selected observable. One thus ends up defining twisted 
fermionic or bosonic fields, the Fourier modes of which are subjected to the following 
(anti) commutation relations: 

SkiAiO-kaAa = ^e 2iklAk2 ak 2 A 2 5-kiAi , (6) 

alxAx 2 = m 2ik ^l 2X X^ ( 7 ) 

CaA^ = ^aXa: " (2vr) 3 5(k 2 - kO^Je 2 * 1 ^, (8) 

where hk represents the spatial components of momentum and the {A}'s represent the 
spin degree of freedom. Incidentally, this structure of twisted (anti)-commutation relation 
can also be obtained in the non-relativistic domain, if one considers the twisted action of 
the Galileo group As has been shown in the paper just mentioned, one has to just 
restore the symmetry corresponding to spatial rotation (SO(3)) here, as the Galileo boost 
generator does not get affected by the twist. Also the wedge symbol (A) inserted between 
ki and k 2 denotes 

ki A k 2 = h u e^k 2j . (9) 

Finally, rj = ±1 as dictated by Boson or Fermion statistics. As we are interested in 
fermionic systems r\ — — 1. These operators go over to usual Bose/Fermi ladder operators 
in the limit 9 — > 0: a — > a. In fact, it has been shown in [5] that a and a are related as 

6 = e (fM^,)a. 



3 A brief review of twisted (anti)-commutation rela- 
tion 

It has been shown earlier that in a non-relativistic system, the presence of spatial non- 
commutativity given by # tJ spoils the SO (3) rotational symmetry. This can be seen easily 
from the fact that the vector = {6i}, dual to the second rank anti-symmetric tensor 9^ 

Bi = \e ijk 9 jk , (10) 
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fixes a direction in the 3-dimensional space. Thus the only surviving symmetry is the 
SO (2) rotation around the direction. However, applying the Drinfeld twist to the Hopf 
algebra, the entire SO (3) symmetry can be restored. This is done by considering a twisted 
co-product, Afl(My), corresponding to SO(3)-generators My, with action defined as 

A e (M ij ) = Ff 1 A (Mv)Fo, (11) 

where 

A {My ) = My (g) 1 + 1 (g) My (12) 

is the co-product in the commutative case (9 = 0) and F e = e~^ e,3Pl ® Pj is the twist 
operator. The corresponding co-product for the rotation group SO (3) in the commutative 
and non-commutative cases are 

A Q (g)=g(g>g, A e (g) = F e - 1 A (g)F e . (13) 

Note that the star-product, defined in (j2J), can now be written alternatively as 

(a */3) (x) = m e (a (g)(3) = m (F e a (g) (3) (14) 

where mo represents the composition of a pair of functions in the commutative (9 = 0) 
case, which is nothing but point- wise multiplication of these functions,i.e, a(x)(3{x). 

It can be easily seen at this stage that the usual projection operator for a two-particle 
system, P = |(1 ± r ), involving the flip map t (o:®/3) = /3(g) a, and which projects 
onto the symmetric (anti-symmetric) sub-space describing bosonic (fermionic) statistics, 
does no longer commute with Ag: [Ag, P ] ^ 0. This implies that one must twist the flip 
map as 

Te = F e - 1 T F e = F e - 2 T 0l (15) 

so that the corresponding projection operator Pg = |(1 ±tq) = Fg 1 PqFq commutes with 
Ag 

[A e ,P e ] = 0. (16) 

This ensures that the new statistics, i.e. the flip operator Tg, remains super-selected, i.e. 
it commutes with all Galilean generators, and defines twisted bosons or fermions. We 
now provide a heuristic argument how the algebra involving the corresponding creation 
and anihilation operators also get modified. 

To begin with, let us apply the twisted projection operator Pg on the tensor product 
of two momentum eigen- states \k > \ l > and make use of f|T5l) : 

Pe(\k >(g) \l >) = \{\k > \l > +r]Fg- 2 \l ><g\k >). (17) 

Up to an overall phase, this can be re-written in a symmetric or anti-symmetric form as 
in the 9 = case 

Pe(\k >(g) \l >) = j khl [\{\k,l » k >>)], 

\k,l »= e~ ikAl \k >(g)|/ > . (18) 

Now identifying a fc a ; |0 >= Pg(\k > \l >) it easily follows that a k a l = T]e 2%kAl aia k . The 

other phase e~ 2lkAl , occuring in a k ai of eq. (jHJ), can be easily understood from the fact 
that the anihilation operator a\ is associated with a momentum (— I), in contrast to the 

operator a], for which the associated momentum is (+/)• 
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4 Non-commutative degenerate electron gas 



To begin with, the Hamiltonian operator for the electronic system in the second-quantized 
formulation can be written as 

H el = J d 3 xft(x) * Tty(x) + 1 J J d 3 xd 3 y^(x) * * V(x, y) * ^(y) * ^(x). (19) 

Before we explain the terms and notations, first note that in order to get a well-defined 
thermodynamic limit, we introduce a box of volume V = L 3 , containing N particles, as 
a regulator so that when V — > oo and iV — > oo, the density (N/V) is held fixed. This 
implies that the field operators 

Vi(x) = E^kA(x)a kA , (20) 

k,A 

occuring in ( fT9l) act on the Fock-space of states obtained by superposing anihilation 
operators in terms of the single particle wavefunctions 

<Mx) = ^=e< kx r/ A (21) 



V 

with k taking the discrete set of values k = with n representing a triplet of positive 
or negative integers, arising from the imposition of periodic boundary conditions in the 

box, and r]\ stands for two spin functions, = ( q ) > ^ = ( 1 ) ' Furthermore, the 

arguments x are now the usual c-numbered coordinates, so that they compose through *- 
products. In this form the action is hermitian, as follows from the fact that (f*gy = g^*f^ 
[7]. Here we would like to clarify some conceptual issues. These are related to the 
concept of a box in NC space. First of all, at the level of operators x^ p satisfying JT]) the 
concept of a well delineated boundary does not make sense, as the uncertainty relation 
| Ax* 1 1 Ax- 7 1 > following from ([I]) prevents this. However, at the level of commuting 
x l 's, when the Moyal star-product is being used, one can presumably introduce a boundary 
in the Moyal plane if the "location" is interpreted in the sense of "expectation value". 
Secondly, for particles confined within a box, a typical single-particle wave function has 
a support only within the box, where for a momentum eigenstate this can be taken to 
be in the plane- wave form (T2~T!) . However, these wave functions with a compact support 
do not belong to the Schwartz Class, as mentioned earlier. The best that we can do at 
this stage is to extend the domain of definition of the function from the box to the entire 
three-dimensional space through periodicity and then compose any pair of such plane- 
wave functions through the star-product. The resultant function, which incidentally is 
also of the same form (T2~Tj) up to a phase factor, is then restricted within the box. One 
can expect that this prescription will yield the desired result in the thermodynamic limit 
L — > oo. 

The kinetic and interaction energy operators are denoted by T(x) and V(x, y), re- 
spectively: 

p 2 -|ii|x-x| 

V{^y) = ~ r . (23) 
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This has to be augmented by the Hamiltonian of the positive inert background having 
the particle density n(x) 

H b = — drxdTx' w , ' — ; (24) 

2 J |x — x'| 

and the Hamiltonian 

N 



tlel-b — — 6 2^ d X i i ' V^" 3 / 

i=r l x ~ x *l 



representing the energy between the electrons and positive background. With this the 
total Hamiltonian for the system becomes 

H = H el + H b + H el _ b . (26) 

Here we have introduced a screened Coulomb potential through an exponentially 
damping factor having a regulator fi of dimension [L] -1 . This renders the integrals ap- 
pearing in H b and H e i_ b finite in the thermodynamic limit. It should, however, be kept in 
mind that the limit \x — > should be taken after the L — > oo limit, so that one can ensure 
fi^ 1 < L at each step of the computation. This also facilitates the shifting of origin of 
integration, as dictated by convenience, as the surface terms are negligibly small in this 
limit. 

Taking the distribution to be uniform n(x) = N/V, we can now compute the pair of 
non-dynamical (inert) terms 

*=2 e V?' (2T) 



H d - b = -e^% (28) 



V 

where we have made use of the translational invariance. As expected, these are c-number 
terms. The total Hamiltonian (T26]) thus reduces to 

1 9 N 2 4ir . , 

= "2 e V? + H " (29) 

with all the interesting physical effects being buried in H e i. 

We therefore turn our attention to H e \. To begin with, let us consider the kinetic term 
in f[T§|) first. This can be simplified as, 



c/ 3 x^ f (x) * 7ty(x) = E ^ — a k , A ,a kA / rf 3 x^/ V (x) * VkA(x). (30) 

kAk'A' ^ m 



Using the definition of the *-product given in (jSJ), we can easily see that this brings in an 
exponential factor e~* k Ak , as the ^-operator is sandwitched between a pair of plane-wave 
states of the form f[2~Tl) . However, the integration yields ^kk'^AA', thereby reducing this 
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^-dependent exponential factor to identity. As far as this kinetic term is concerned, this 
therefore yields the same result as in the commutative (8 = 0) case: 



r ~ „ ~ h 2 k 2 

/ c/ 3 X^ f (x) * Tlp(-X.) = T, «kA fi kA, (31) 

kA Zm 



which can be interpreted as the kinetic energy of each mode multiplied by the correspond- 
ing number operator. Almost the same thing also happens for the potential energy term 
as all ^-dependent exponential factors coming from the Moyal star-product reduces to 
the identity here as well. Seemingly, the non- commutative structure of space, encoded 
in the Moyal star-product, plays no role, at least to first order in perturbation theory. 
This was also found in [3]. There the effects of non-commutativity only showed up in 
second order, which is not unexpected as the exchange correlations of a non-commutative 
theory will generacilly be different from those of a commutative theory. However, as we 
proceed to show now, a ^-dependence, stemming from the twisted anti-commutation rela- 
tion, survives to give a NC modification to the ground-state energy already to first order 
in perturbation theory. 

To this end consider the potential term in ffl9|) . First of all this requires a proper 
interpretation. Given that x and y represents the position co-ordinates of two distinct 
particles in a pair, we must have y v ]* = 0. Consequently, only functions involving 
coordinates of a single-particle will compose through the *-product. This implies that the 
appropriate way of writing this potential term is 

d 3 a^t( x ) * ( / d 3 y^(y) * V(x, y) * ^(y)) * ^(x), (32) 



where the y-integral is performed first, after composing the three functions of y, with x 
being held fixed, to yield a function of x. With this interpretation, we can write 

/ / d 3 xd 3 yij\x) * ij\y) * V(x,y) * ^(y) * ^(x) = - J2 5 k 1 A 1 «k 2 A 2 ak 3 A 3 ak4A4 x 

WW 

< k 1 A 1 k 2 A 2 |y|k 3 A3k4A4 >, (33) 

where the matrix element is given by 

2 

< kxAxksAal F|k 3 A3k 4 A4 >= ^ / d 3 xe- l ^- x VXl (iy * ( J d 3 ye~^ Vx2 (2)U 

f^4*e 4k ^ A4 (2)) * e^- x 77 A3(1) . (34) 

l x - y| 

Making use of the identity 6 i~~r = t^W / J^+k 2 e^*^ 0we can simplify this matrix 
element to get 

< k 1 A 1 k 2 A 2 | F|k 3 A3k 4 A4 >= J ^^e- jklAk3 e- ik2Ak4 5 AlA3 5 A2A4 x 

0(kl+k 2 ),(k3+k 4 )) (35) 



1 This can easily be seen to follow by noting that the screened Coulomb potential is a Green's function 
of the Laplacian augmented by a mass term: (—V 2 + /i 2 )( e a ■ ) = 47r<5 3 (x) 
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where k gets restricted to k = k x — k 3 = k 4 — k 2 , so that k, occuring in the Fourier trans- 
form of the Screened Coulomb potential, can be identified with the momentum transfer. 
On the other hand, the Kronecker delta involving momenta levels enforces momentum 
conservation. We therefore convert integration over to a sum by the standard re- 
placement / d 3 k — > (^) 3 so that the parity with other momentum variables can be 
restored. We thus get, 

2ttp 2 1 

< k 1 A 1 k 2 A 2 |F|k 3 A 3 k 4 A4 >= -^E^-TI^ tklAk3e " k2Ak45 A 1 A3 5A 2 A 4 5(k 1+ k 2 ),(k 3+ k 4 ).(36) 

k ft 

The potential energy operator becomes 
2 7I -g 2 1 

~TT E 2 i 72 e iklAk3e_lk2Ak4 ^AiA3^A 2 A4^(k 1 +k2),(k3+k4)«kiAi«k2A2 5 'k4A4«k3A3, (37) 

V W{ x}ft + k 

where {k} represents the set (k, ki, k 2 , k 3 , k 4 ) on which the summation has to be per- 
formed. The total Hamiltonian can now be written as 

1 oN 2 47r h 2 k 2 „± _ 2ne 2 ^ 1 , k AW , w AW . 

H =-~Z V— +El7-4A«kA + - f7 - E -Tm e 1 3e 2 4J A 1 A 3 5a 2 A4 X 

2 V fx 2 2m ^ { kH A} ^ 2 + ^ 2 

^(k 1 +k 2 ),(k3+k4) a k 1 Ai a k2A2^k4A 4 ak3A3- (38) 

The electrical neutrality of the system makes it possible to eliminate \i from the 
Hamiltonian. The last term of the Eq. ([58]) can now be re-cast using the momentum 
transfer k as 

Ittp 2 1 

-FT" E -^e-^e^^4 +k3Ai 4 4 _ kA2 S k4A2 S k3Al , (39) 

V [kk 3 k 4 A 1 A 2 ] ft +K 

where the two spin summation have been evaluated with the Kronecker deltas. At this 
stage, it is convenient to seperate the above expression into two terms, refering to k ^ 
and k = 0, respectively, 

9-7TP 2 ' 1 

Z7rC \^ 1 p -ikAk 3g ikAk 4 ~t ~t ~ ~ , 

y 2 I jU2 e e a k+k 3 ,Ai a k 4 -k,A2 a k4A2 a k 3 A 1 + 

[kk 3 k 4 AiA 2 ] " 

27re 2 ^ 1 _+ _+ 

— T7 - 2^ — a k 3 Ai a k4A2 a k4A2 a k 3 A 1 (40) 

K [k3k4A!A 2 ] ^ 

where the prime on the first summation means the k = term is omitted. The second 
term may be rewritten with the discrete version of twisted anticommutation relation given 
by Q), ([8]). In this discrete version, (Q remains unaffected, while one has to just replace 
(27r) 3 5 3 (ki — k 2 ) — > 5kik 2 i n &■ Using this the second term takes the form, 

^ 7Ce ' 2 1 If-t = X at X r ~ \_„2ik 3 Ak 4 „2ik 4 Ak 3 
2^ T^-l^sA^kaAi^Aa^Aa - OkaA^kaki^AxAa^AaJ^e e 



^ [kakiAxAa] ft ^ 



— ^(iV 2 -iV), (41) 
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where 



N — ( d 3 xn(x) = E a% = E a\a r = E™r = j d 3 x^(x)tp(x) (42) 

«* |V» IV" (*» W 



represents the number operator. Here too the effect of non-commutativity disappears. 
H Since we always deal with states of fixed N, the operator N may be replaced by its 
eigenvalue N, thereby yielding a c-number contribution to the Hamiltonian 

N 2 e 2 27r Ne 2 2n . . 

— 7? - ~7- < 43) 

The first term of the above expression cancels the first term of the Hamiltonian in ( 1291) . 
The second term represents an energy —2ne 2 (Vfi 2 )~ 1 per particle and vanishes in the 
proper physical limit: first L — > 00 and /i as discussed earlier. Thus the explicit /i~ 2 
divergence cancells identically in the thermodynamic limit, which reflects the electrical 
neutrality of the total system; furthermore, it is now permissible to set \x = in the first 
term of (1391) . since the resulting expression is well defined. We therefore obtain the final 
Hamiltonian for a bulk electron gas in a uniform positive background 

h 2 k 2 2ire 2 ' 1 

^ = £"2^4*^ + — £ £ pe- lkAk3 e lkAk4 4 +k3iAl 4 4 _ k)A2 a k4A2 a k3Al . (44) 

kA kk3k4 A1A2 

We can now introduce the inter-particle spacing r , defined through V = ^nr^N and 

measure it in units of the Bohr radius ao = -^-2, thus yielding a dimensionless variable 

10 

With r as the unit of length we define the following quantities 

V = Tq V, k = r k, k 3 = r k 3 , k 4 = r k 4 (45) 
and thus we obtain the following dimensionless form of the Hamiltonian operator 

H = -^(£ ^ 2 4 A a £A + ^ E £ p4 + k 3 ,A 1 4 4 -k, 2 ^A 2 ^A 1 e--o kA ( k 3- k 4 )).(46) 
a ° r s k A 1 v [kk 3 k 4 ] A1A2 

In the limit r s —>■ 0, corresponding to the high density limit (tq —>■ 0), the potential energy 
becomes a small perturbation even though it is neither weak nor short ranged, while the 
leading term comes from the kinetic-energy term. Thus the leading term in the interaction 
energy of a high-density electron gas can be obtained with first order perturbation theory. 
In the high-density limit, we can therefore separate the original dimensional form of the 
Hamiltonian (1441) into two parts: 

^0 = £ ^aUa, (47) 

kA 



Hi = ^T £ E^U.A^L-k.A^k^akaA.e-^-^, (48) 

[kk 3 k 4 ] AiA 2 

2 As has been mentioned earlier it was shown in reference [B] that the twisted and untwisted ladder 
operators are related by a unitary transformation, thereby preserving the number operators 
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where Hq is the unperturbed Hamiltonian, representing a noninteracting Fermi system, 
and H% is the (small) perturbation. Correspondingly, the ground state energy E may 

be written as E° + E 1 + , where E° is the ground-state energy of a free Fermi gas, 

while E 1 is the first-order energy shift. Since the Pauli exclusion principle allows only 
two fermions in each momentum eigenstate (hold even in noncommutative space jl, one 
with spin up and one with spin down, the normalized ground-state \F > is obtained by 
filling the momentum states up to a maximum value, the Fermi momentum pp = Tikp. In 
the thermodynamic limit one can again replace summation by integration, so that kp can 
be determined by computing the expectation value of the number operator in the ground 
state \F > 

N =< F\N\F >=J2< F \^\ F >= E & ( k F -k) = (37r 2 )-V4 = N (49) 

kA kA 

where Q(x) denotes the step function 

Q(x) = 1 for x > 0, 0(x) = for x < 0. (50) 
Equivalently k F can be expressed as 

, ,3n 2 N i 971. i _! 

k F = (— — ) 3 = (^-) 3 r \ (51) 
Now the expectation value of Hq may be evaluated as, 

=< F\H \F >=^Y,k 2 < F\n kX \F >= ^ £ fc 2 6(A; F - fc) 

kA kA 

= y 7 — ^ / d 3 kk 2 Q(k F -k) — ^iV = t^(^) 3 - 52 

In a free Fermi gas, the ground-state energy per particle E^/N is | of the Fermi energy 
e F = h 2 k 2 F /2m. We now compute the shift in the ground-state energy to first order in 
perturbation theory: 

E {1) = <F\H 1 \F> 

?TTP 2 ' 1 

= "IT ^ ETj<f|4 +k3 ,A 1 4 4 - M2 «k 4 A 2 a k3Al e- !k ^- k ^>. (53) 

kk 3 k 4 AiA 2 

For a non zero matrix element, all the states |k 4 A 2 >, |k 3 Ai >, |k + k 3 , Ai >, |k 4 — k, A 2 > 
should be occupied in \F >. Also, since k = is excluded, we must have the pairing 

k + k 3 , Ai = k 4 , A 2 



As one can see that the Pauli's exclusion principle remains valid in the momentum space, as can 

be seen from (6) that in the case of r\ — * — 1, a k a k = 0. Consequently, the concepts like Fermi label 
etc. remains valid in momentum space even in presence of non-commutativity. However, this does 
not remain true in configuration space, as ip^ (x)ifj^ (x) does not vanish necessarily as follows from the 
relation ip(x)ip(y) = J d 2 x'd 2 y'Tg(x, y, x', y')4>(y')ip(x')] 9 ^ 0[4]. Indeed it has been shown there, that 
repulsive statistical potential between a pair of identical (twisted) fermions can saturate to a finite value 
at coincident points, thereby violating Pauli principle in real space. 
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so that the matrix element becomes 



<Wk 3 ,k4^AiA2 < -^l a k+k3,Ai a k4-k,A 2 a k4A2 a k 3 Ai|-^ > • 

The exponential factor in the previous expression for (pD again becomes the identity 
as k + k 3 = k 4 (z.e. k = k^ — k 5 ). Finally, using the twisted commutation relation given 
by ([8]) and also considering the fact that the term k = is excluded from the sum, 
the above expression becomes 

- 5u +ksM 5 XlX2 < F|n k+k3iAl n k3Al e- 2ik3A(k+k3) |F > 
= -<Wk 3 ,kAiA 2 e(A; F - |k + k 3 \)Q(k F - fc 3 )e- 2lk3Ak . (54) 



From this follows 



9-7TP 2 1 

E^ = --^- £ ^Q(k F -\k + k 3 \)Q(k F -k 3 )e- 2 ^ 

V [kksAx] fc 



E(^) 2 / d 3 kd%-Q(k F - |k + k 3 |)0(A; F - fc 3 )e- 



Note that we have again included the k = term at the level of integration, as being 
a set of measure zero, it does not contribute anything to the integral. Now taking the 
summation over Ai the above expression becomes 

= -^p- J d 3 kd 3 k 3 k~ 2 e(k F - |k + k 3 \)B(k F - A; 3 )e- 2ik3Ak 

It is convenient to change variables: k 3 — > k 3 + |, which reduces the above expression 
into the symmetrical form 

= / ^k 3 Q(k F - |k 3 + Ik|)6(^ - |k 3 - ik|)e-^%. (55) 

The computation of this integral can now be made simpler by orienting the 3rd axis along 
the dual 6 from eq.(10), so that only 9 12 ^ (see appendix A). Although a bit lengthy 
the integral can then be computed in a straightforward manner (see appendix B) to yield 

Eil) _ _ levrw ~ {-mm^z (56) 

(2tt)« ^^(i + r) (i + 2r)(r(l + r))»- W 

This is a infinite series in 9. Now as 9 is of the order of Plank length-scale, hence very 
small, we sum the series up to the first contributing term in 9, neglecting all other higher 
order terms. This gives 

(2tt) 6 m 54 ; 2a LV 4 ; 2vrr s v 256 ; ( a 2 J r 5 / V ; 

Thus the ground-state energy per particle in the high density limit is given approximately 
by 

E e 2 r 2.21 0.916 , 9 N2 0.2302 

+ (-2)—^ + 1- ( 58 ) 



N r a ->o 2a 1 r 2 s r s a 2 , r 5 s 
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As can be seen from the above expression, the ground-state energy has 9 corrections and 
deviates from the commutative result. The ground-state energy per particle is increased 
by an order of (~ (A) 2 ), which is a dimensionless quantity. Note, however, that this 9 
dependency has its roots in the twisted anti-commutation relations, rather than the non- 
commutative structure of space per se, as the 9 dependency arising from the star product, 
which encodes the non-commutativity of space, always dropped out. This is true to first 
order in perturbation theory. However, as was shown in [3]) both these effects play a role 
to higher order in perturbation thery. 

It is also quite clear from the above expression that when r s becomes very small, 
i.e. in the high-density limit, the effect of spatial non-commutativity on the ground- 
state energy of the degenerate electron gas becomes more significant. Taking the non- 
commutativity parameter (9) to be of the order of plank-length it can be seen that the 
effect of non-commutativity in any terestrial experiments may not be found, but in the case 
of astrophysical objects, where matter density is very high, the effect of non-commutativity 
may be found by experiments. 

5 Conclusions 

We have shown that, in contrast to [3], the ground-state energy of a non-commutative 
degenerate electron gas in a neutralising background acquires non-commutative correc- 
tions to first order in perturbation theory. These corrections arise from the modified two 
particle correlations resulting from the twisted anti-commutation relations. All 9 depen- 
dency arising directly from the star product dropped out to this order in perturbation 
theory. Our final observation is that any observable effect, as far as the energy shift is 
concerned, can only be obtained when the system is extremely dense - a situation that 
can presumably arise only in an astrophysical setting. 
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7 Appendix A 

Here we would like to demonstrate a problem in 3-dimension with spatial non commuta- 
tivity can be reduced effectively to a problem of planar non commutativity by making an 
appropriate SO (3) transformation. 

To that end, let us first recall that the vector = {9i} dual to the non-commutative 
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antisymmetric parameter 9 %3 is defined as in eq.(10), 

1 







Now this 0-vector will be pointed in an arbitrary direction. So we can parametrise it as 

/ sin a cos (3 \ ( 9 23 \ 
6 = 9 sinasin/? = 9 31 . (59) 
^ cos a / \9 12 J 

We would now like to rotate the coordinate axes XYZ to X'Y'Z' by making an SO (3) 
transformation in such a manner that Z' axis now becomes parallel to the 6 direction. A 
convenient choice for the element of SO (3), doing this job, is obtained by taking products 
of two successive SO (2) rotation matrices around Y and Z-axes as 



R = 



I cos (3 — sin (3 
sin (3 cos (3 

V o oi 




cos a cos (3 — sin (3 sin a cos (3 
= | sin (3 cos (3 cos (3 sin a sin (3 
— sin a cos a 



e SO (3) 



(60) 



As can be easily checked that the action of the i?-matrix on the canonical orthonormal 
basis, 




;e 3 




(61) 



associated with the original XFZ-frame, yields following the orthonormal basis e- = Ref. 

(62) 



/ cos a cos (3 \ 




(-sin(3\ 




( sin a cos /3 \ 


cos a sin (3 




cos (3 




sin a sin /3 


\ — sin a ) 




\ 1 




V cos CK / 



for the X'y'Z'-frame, as can read-off easily from the respective columns of the i?-matrix 
(). Particularly, one can verify that e' 3 is indeed proportional to the 6- vector. 

Now the position vector r can be expanded in either of these basis as, r = x*e, = 
x' 3 e'j. Making use of the fact that e- T e^- = 5ij, we can obtain x' % in terms of x 3 as 

x n = x 3 \e lT R~ l ej). Writing elaborately, this yields 



fx' 



S] I ■ » \ / X \ 

y 



I cos a cos (3 cos a sin (3 

— sin (3 cos a 
\ sin a cos (3 sin a sin (3 cos a / 

/ x cos a cos (3 + y cos a sin /3 — z sin a \ 
= —x sin (3 + y cos /3 

\ x sin a cos /3 + y sin a sin /3 + z cos a / 

Now, if the original x^s are operator- valued and satisfy 

[x\ x i] =iQ i3 \ 



(63) 
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then one can easily check that 



[x',y']=ie, [y',z'} = [z',x']=0. (64) 

Thus the problem is essentially reduced to planner non-commutativity. This result is 
some-what expected, as the matrix {%} is degenerate, as it is an odd dimensional (3 x 3) 
anti-symmetric matrix. And we have achived this by orienting the z'-axis of the rotated 
coordinate axes in the direction of 6- vector introduced in (10). 



8 Appendix B 

Here we provide some of the basic steps leading to the evaluation of the integral 

r d 3 k 1 1 

1 = J -j^d%Q(k F - |kg + ^\)Q(k F - |kg - -k\)e^ k ^ k >. (65) 

These functions indicate that the exponential function has to be integrated first over 
the overlapping region 'i?' of the two spheres of equal radii k F centred at ±(|k) with 
respect to the variable k3 holding k fixed. The resultant function of k is then integrated 
over the solid sphere |k| < 2k F . Considering the region R + (k 3z > 0) initially, we can 
write the corresponding integral over k 3 with proper limit as, 



L+= f d 3 k 3 e- ik ' Mal where a* = 9 ij kj 

J Rj- 



dze~ iaz / * dye~ iay / dxe 

Here we have substituted k 3x = x, k 3y = y and k 3z = z for convenience. 
Upon simplification, this takes the following form, 



L+ = ]_ [ {kF J dze -io?z f +p dye - ia S sin ( a i £?Zy2) where p = J kl - (z + -) 2 . 
cr Jo J-p v V 2 

Substituting the above expression for L + in (62), one gets for the corresponding I + 

I+ = J 1^ L+ = 1 Si F J dz Jo d y Cos ( a2 y) Sin (^\/p 2 - y 2 )l 

where we have taken 9 13 = 8 23 = 0, by orienting the 3-axis in k frame in the direction of 
0, where 9{ = \tij\$jk (see (10) and Appendix A) so that the only surviving componenets 
of 6 l i are 9 12 = —9 21 = 9 consequently a 3 = and a 1 = 9k 2 and a 2 = —9ki. The above 
integral can now be expressed in terms of Bessel functions [9] 



as 



2^=0 2 2r r\T(v + r + 1) 
f d 3 k /"( fc f-|) , p . r. — . — 
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Now we can integrate term by term by expanding series of Bessel functions. By making 
use of the integral 

/•(**■- k r+1 r+1 <7f— 1 V h 

/ dz{ki - (z + yy +i = (k F r+v m ii (i _ (-f 

•<o 2 ^ 2j + 1 2fcp 



one gets after a lengthy but straight forward computation 

~ (_l)r( fc 2 g)2r r+1 r+lQ.^j r rQ^f 

+ " F h rKr + l)!(2r + 1)! ^ (r + J + 1) (2/ + 1) ' 



The last two factors involving series summation can now be carried out exactly to get 



^(2/+l) 2r(f+r)' 
£S F+7+1) " r(f + r) 



With this J + takes the following form 



/+ = 4 " E n 2fl + r)fl + 2r)frf| + r^- (66) 



r=0 



A similar result is obtained for J_, when one considers the region i?_ corresponding to 
k 3z < 0. Consequently one has I = 2I + = 2J_. Substituting this back in (54), one gets 
the desired expression of E®. 
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